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ABSTRACT 

Many hydrogen deficient stars are characterised by surface abundance patterns that are hard to reconcile with 
conventional stellar evolution. Instead, it has been suggested that they may represent the result of a merger 
episode between a helium and a carbon-oxygen white dwarf. In this Letter, we present a nucleosynthesis 
study of the merger of a 0.4 Mq helium white dwarf with a 0.8 Mq carbon-oxygen white dwarf, by coupling 
the thermodynamic history of Smoothed Particle Hydrodynamics particles with a post-processing code. The 
resulting chemical abundance pattern, particularly for oxygen and fluorine, is in qualitative agreement with the 
observed abundances in R Coronae Borealis stars. 

Subject headings: Nuclear reactions, nucleosynthesis, abundances — stars: AGB and post-AGB — stars: white 
dwarfs 



1. INTRODUCTION 

Hydrogen deficient stars with high carbon abundances 
(with the exceptions of the PG 1159 and CSPN WR stars) 
can be divided into three main categories according to their 
effective temperatures: Hydrogen Deficient Carbon (HdC) 
stars, R Coronae Borealis (RCB), and Extreme Helium (EHe) 
stars. They are often considered to be at different evolution- 
ary stages following a common origin, as evidenced by their 
similar atmospheric abundances (Pandey et al. 2004). 

The origin of these stars is of particular interest because 
of their unique chemical composition. As well as being en- 
riched in carbon and oxygen relative to solar abundances, 
they also exhibit enrichment in nitrogen, neon, fluorine, and 
some s-processed material. These abundance patterns are 
hard to reconcile with their expected initial composition, re- 
quiring alternative scenarios for their origin. Two leading 
theories are the Final Flash (FF) and the Double-Degenerate 
(DD) scenarios. The FF scenario attributes the formation to a 
late He-shell flash in a post-Asymptotic Giant Branch (post- 
AGB) star as it evolves towards the white dwarf state. The 
flash would move the star back onto the AGB part of the 
Hertzsprung Russell diagram, while mixing s-process ed ma- 
terial and He- b urning ashes to the surfa ce (Fuiimoto 1 19771 : 
iRenzinil 119791: IWerner & Herwig 120061) . In the DD sce- 
nario, the merger of two white dwarfs (He-i-CO), with a to- 
tal mass below the Chandrasekhar limit, through disruption 
and accretion of the lighter star, gives rise to a single object 
whose composition is an admixture o f the two white dwarfs 
(IWebbinklll984l: lAsplund et al.l |2000'). This dense object is 
expe cted to be surrounded by a hot corona and an accretion 
disk dLoren-Aguilar et aTl 120091) . Surface abundance analy- 
ses of hydrogen-deficient s tars seem to favour the latter sce- 
nario dSaio & JefFervll2002l) . although the possibility of some 
objects being formed through the FF scenario cannot be ex- 
cluded (Pandev & Lambert 20 11). 

When considering the formation of hydrogen deficient ob- 
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jects through the DD scenario, the following question can 
be raised: does nuclear processing occur during the merg- 
ing event? If no nuclear burning occurs — the "cold" merger 
hypothesis — the observed abundances might be reproduced 
through the partial mixing of material previously processed 
during the AGB phase of the CO white dwarf progenitor 
However, recent calculations based on this scenario failed to 
quahtatively match the overall abundance pattern observed 
(Jeffer v et al.l 1201 Ih . Another possibility involves a "hot" 
merger event, in which signi ficant nuclear processing occurs. 
Indeed, IClavton et"al] (120071) suggested that this scenario may 
account for the high abundances of '^O observed in some 
stars. Unfortunately, no detailed nucleosynthesis study dur- 
ing white dwarf mergers has been published to date. 

In this Letter, we present a detailed study of the nucleosyn- 
thesis accompanying white dwarf mergers. Calculations are 
based on post-processing analyses using temperature and den- 
sity versus time profiles obtained from Smoothed Particle Hy- 
drodynamics simulations of the merger event. The structure 
of the manuscript is as follows: in section |2] we summarise 
the observed abundances in hydrogen deficient objects. The 
details of previous nucleosynthesis studies and of the calcu- 
lations presented in this work are given in sections [3] and HI 
respectively. Results are discussed in section |5] and the main 
conclusions are reached in section|6] 

2. ABUNDANCES IN H DEFICIENT STARS 

All current observed surface abundances of hydrogen de- 
ficient stars have been collected recently by I Jeff'ery et alJ 
(1201 llrl to which the reader is referred for details. Here we 
summarise the main observational features reported in that 
work. 

Carbon is enriched in all observed hydrogen deficient stars, 
irrespective of their iron content, with overabundances rela- 
tive to solar of the order of 0.7 dex. Additionally, the inferred 
'^C/'^'C ratios are very large (> 500), indicating helium burn- 

' with the exception of '-C/'^C and '^0/'**0 ratios, and abundances of 
Sc, Mn, Co, Cu, and Sr 



ing in the star. Nitrogen is also enriched in most EHe and 
RCB stars. However, in this case it is correlated to the iron 
abundance, suggesting that nitrogen is enriched only by CNO 
cycling prior to the merging event, as can be shown by sum- 
ming the initial carbon, oxygen, and nitrogen abundances of 
the star. Oxygen poses an interesting problem since its abun- 
dance exhibits a large spread. Regardless of this fact, oxygen 
appears to be both enriched and independent of metallicity. To 
account for the overabundance of oxygen a dredge-up event, 
by which oxygen is injected from the CO white dwarf, seems 
to be required. Because of the difficulties in distinguishing be- 
tween isotopes in hotter stars, oxygen could also be produced 
in the form of '^O, which has been found to be highly enriched 
in a few stars dClavton et al.ll2007l) . Fluorine, in turn, is highly 
enriched by 2^ dex in all the observed stars. This enrichment 
could possibly arise from processing during the AGE phase of 
the CO white dwarf progenitor through a rather complicated 
process in the helium inter-shell zone (Lugaro et al. 2004). Fi- 
nally, s-process elements, particularly yttrium and zirconium, 
have also been obs erved to be enriched by up to a factor of 50 
in some EHe stars iPandev et al.l (120041) . The source of these 
elements is probably dredge-up of helium inter-shell products 
of the AGB star that subsequently evolved into the CO white 
dwarf. While several models can match some of these ob- 
servational features to an order of magnitude level, the cur- 
rent challenge is to reproduce the overall abundance pattern 
simultaneously with a detailed model of the merger event. 

3. PREVIOUS NUCLEOSYNTHESIS STUDIES 

Most abun dance studies of white dwarf mergers per- 
formed so f ar (Saio & Jefrer3|200l|^deiL& Lamberl|20ll 
iJeiferv efaJ. 2011) consider the cold merger paradigm, in 
which no nucleosynth e sis oc curs. In particular, the recent 
work by iJefferv et al.l ( 12011b started with the abundances 
present in the white dwarfs as resulting from previous evo- 
lution, particularly the AGB star evolution of the CO white 
dwarf progenitor. Calculations also included helium and hy- 
drogen buffers on the surfaces of the white dwarfs, remaining 
from previous evolution, which contained material that had 
been heavily processed. By using a recipe for mixing the ma- 
terial in the white dwarfs during the merger, they were able 
to account for most of the abundances (or at least the trends) 
observed in EHe and RCB stars. The most notable exception, 
however, was the abundance of fluorine. Although fluorine 
can exist in the helium buffer of the CO white dwarf, it was 
found not to be abundant enough to account for the very high 
observed values. 

Hot merger nucleosynthesis calculations were performed 
by IcTavton et al. (2007) upon the discovery of very low 
i6q^i8q j-^fJQg jjj some HdC and RCB stars. These studies 
relied on a one-zone parametric model to simulate the con- 
ditions needed to burn helium during a merger event. Their 
simple model was able to account for the enrichment of '^O 
in these environments. When an admixture of hydrogen was 
added, however, the '**0 was rapidly destroyed through proton 
captures while "'O was created through the '^CCa, n)'^0 re- 
action, thus restoring the '^O/'^O ratio closer to solar values. 
Hot merger nucleosynthesis was also studied by IJefferv et alj 
(12011) in the absence of mixing of material from the CO 
white dwarf. They were able to account for overabundances 
in carbon a nd oxygen, but limitations in the SPH calculations 
provided bv lLoren-Aguilar et alj (l2009h prevented them from 
making reliable predictions for other nuclei. Detailed nucle- 
osynthesis studies using reaUstic hydrodynamic models are 



therefore required to fully understand the production of ele- 
ments in these environments. 

4. INITIAL MODELS AND INPUT PHYSICS 

The models used as input for the nucleosynthesis calcula- 
tions reported in this paper h ave been presented elsew here 
dGuerrero et al.ll2004t iLoren-Aguilar et al.ll2009l |20T0). and 
rely on a Smoothed Particle Hydrodynamics (SPH) La- 
grangian simulation to follow the evolution of two coalesc- 
ing white dwarfs. We refer the reader to these works for the 
details of the dynamical behaviour of the merger process, but 
we describe here the structure of the final remnant, which is 
relevant for the calculations reported below. In essence, the 
final remnant of the merger process consists of a central hot 
white dwarf which contains most of the mass of the primary. 
On top of it a hot corona can be found, which contains most 
of the mass of the disrupted secondary and a small admixture 
of the primary white dwarf. Finally, surrounding the com- 
pact remnant, a rapidly rotating Keplerian disk is formed. It 
is also worth noting that little mass is ejected from the sys- 
tem. We emphasise that, for the sake of conciseness, we only 
describe the results obtained when the masses of the merging 
white dwarfs are 0.4 and 0.8 Mq with SPH particle masses of 
2.6 X \Q-^Mq and 5.3 x IO^^Mq, respectively. 

A network containing 14 nuclei was employed in the SPH 
simulations to account for the energy released by nuclear re- 
actions during the merger (mostly a-capture reactions), while 
the detailed nucleosynthesis reported here is calculated for 
tracer particles using a 327 nucleus post-processing network 
ranging from hydrogen to gallium. Energy generation rates 
between the two networks differs by less than one percent for 
our models. The thermonuclear reaction rates are adopted 
from the REACLIB database (Cyburtetal. 2010) with up- 
dates on experimental rates from I liadis et alj (|2010). SPH 
particles that reside in the hot corona of the final object will 
characterise the "atmosphere" that is observed in a hydrogen 
deficient star The tracer particles used to follow the ther- 
modynamic history of the merger are therefore picked from 
0.005 < RIRq < 0.05 as shown in figure [T] This range in- 
cludes particles that represent everything from the surface of 
the central dense object to the inner edge of the accretion disk. 
A total of 10 000 tracer particles have been used here to ensure 
a representative sample from the two white dwarfs. 

Our simulations take into account the existence of thin 
helium and hydrogen buffer shells, which are predicted to 
survive prior evolution through the common envelope stage. 
Therefore, three different regions (carbon/oxygen-, helium-, 
and hydrogen-rich) are defined in the CO white dwarf, while 
two (helium- and hydrogen-rich) are adopted for the He white 
dwarf. These regions are defined a ccording to the mass con- 
tained. Adopting the notation of ISaio & Jefferyl (120021) — 
who, for instance, referred to the mass of the He shell in 
a CO white dwarf as A^fje'CO — ^^^ different shells are 



defined as follows: M, 



CO:CO 



0.78, Ml 



He:CO 



0.019, 



^H:CO - 0.001, MHe:He = 0.399, M^.j^^ = 0.001. The 
chemical composition of each shel l has been taken from de- 
tailed pre-white dwarf models (Ren edo et al.l l2010) of the ap- 
propriate metallicity, with some additional information from 
Heffery et al. (201 1) for the hydrogen buffers. It consists of 7 
species: 'H, '♦He, '^c, i^N, '''O, '^O, and 22Ne. It is worth 
noting that, although we do not expect to reproduce the ob- 
served s-process abundances in hydrogen-deficient stars with 
only these 7 isotopes, they should be enough to model nu- 
cleosynthesis of light and intermediate-mass elements (up to 



neon). Two metalicity scenarios have been considered: (i) so- 
lar metalicity and (ii) a lower metalicity one, with Z - 10"^. 

To account for convective mixing of the material in the hot 
corona, we calculate the mass-averaged abundances of all par- 
ticles that reside within a defined convective zone at the end 
of the hydrodynamic simulation. Although the SPH calcula- 
tions suggest that the entire hot corona will be convective, we 
adjust the depth of the mixing in our analysis to study the ef- 
fect of this assumption on the final surface abundances. Here, 
we consider two cases, (i) "deep" mixing, in which everything 
in the range 0.005 < R/Rq < 0.05 is mixed homogeneously, 
and (ii) "shallow" mixing, in which the range is assumed to 
be smaller, 0.014 < R/Rq < 0.05. 

5. RESULTS AND DISCUSSION 

The mass averaged abundances obtained from the tracer 
particle nucleosynthesis are compared with the observed 
abundances in figure |2] The plotted solid lines correspond 
to the final abundances obtained under the "deep" mixing 
assumption, while surface abundances obtained assuming a 
"shallow" mixing depth are represented by dashed lines. 
These lines are added purely to guide the eye since our cal- 
culations were only performed for two metalicity cases (z - 
0.015 and z = 10"^). In judging the successfulness of our 
models, it is worthwhile to compare each element to observed 
abundances on a one by one basis. 

Carbon depends strongly on the mixing depth assumed in 
the hot corona of the merger product but appears to be con- 
sistently high in our calculations. Nevertheless, the overall 
trend of carbon abundance with respect to metalicity agrees 
with the observational data. The '^C/'^C ratio for solar met- 
alicity is ~ 2 X lO"*", consistent with the predicted nature of 
helium-burned material, but considerably higher than the ob- 
served ratio. This disagreement could arise from the limited 
number of nuclear species adopted in our calculations. 

The oxygen abundance also agrees fairly well with the ob- 
served data, showing a large range of values depending on 
the mixing depth assumed. This dependence could explain 
the observed scatter in oxygen in EHe and RCB stars, where 
the observed oxygen abundance could depend strongly on 
the nature of the individual merging event. This picture is 
not, however, consistent with the relatively low scatter in car- 
bon abundances. With regard to the isotopic abundances, our 
solar metalicity model for deep and shallow mixing yields 
i^O/i^o = 370 and 19, respectively. The "^O has cleai'ly been 
enhanced in the outer regions of the hot corona with respect 
to solar abundances (where '^0/''^0 ~ 400), although not to 
the extent found in some HdC and RCB stars (C layton etal.1 
|2007). 

Nitrogen follows the expected trend indicating that it is en- 
hanced only through CNO cycling in the parent stars. Al- 
though it is slightly reduced by the '^N(q;,7')'^F reaction, it is 
present in such a high concentration that it is barely dimin- 
ished in the short time-scales involved. 



The overabundances of fluorine obtained in our models are 
high, but not high enough to fully agree with the observational 
data. The lack of sensitivity to mixing depth is an indica- 
tion that fluorine is produced homogeneously in the merging 
event, and not in any particular region of the merger. Fluo- 
rine is produced dominantly through the following sequence 
of reactions: ^^N(a,yy^F(p,ay^O(a,yy'^Ne(fi^y'^F. Enrich- 
ment of fluorine in our models is a great success, suggesting 
that hydrogen deficient stars could be the result of a hot white 
dwarf merger event since it is challenging to reconcile those 
high abundances with a cold merger event. 

Finally, the abundance of neon following the merger is con- 
sistent with the metallicity, indicating that it is not processed 
in either the merger event, or during the preceding evolu- 
tion of the parent s tars. Nei ther our models, or the cold 
merger models of Jef Ferv et~an ( 20 11 ) reproduce the observed 
neon abundance in several EHe stars, which contain over- 
abundances up to about 1 dex. Clearly more work must be 
performed to understand the nucleosynthesis of neon in hy- 
drogen deficient stars. 

In summary, our models of hot white dwarf mergers are 
able to reproduce, on an order of magnitude level, the surface 
abundances of light and intermediate-mass nuclei in hydro- 
gen deficient stars. The model's main success is in its ability 
to synthesise high fluorine ab undances, which ca nnot be pro- 
duced in cold merger models (iJeff'erv et al.ll201 ll) . 

6. CONCLUSIONS 

Mergers of helium and carbon-oxygen white dwarfs are 
currently favoured as the events responsible for the produc- 
tion of hydrogen deficient HdC, EHe, and RCB stars. The 
main question facing us is whether these objects are the result 
of a cold merger, in which no nucleosynthesis occurs, or a hot 
merger, in which at least some of the material is processed 
through nucleosynthesis. 

The results presented in this work show that nuclear pro- 
cessing during the merger event can account for some of the 
observed abundances in these stars, particularly for fluorine, 
which cannot be easily synthesised in the cold merger sce- 
nario. There are, however, still open questions regarding the 
expected nucleosynthesis: neon is difficult to produce at the 
observed levels, while the absolute abundance of carbon is too 
high in our models. Such discrepancies may be alleviated by 
the inclusion of a more detailed initial chemical composition 
with many more species, resulting from state-of-the-art mod- 
els of the preceding AGB phase. Inclusion of such detailed 
abundances is currently under way. 
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Fig. 1 . — Radially averaged temperature (solid line) and density (dashed line) profiles of the final He+CO white dwarf merger product as a function of radius 
(in units of 0.1 solar radii). The temperature profile clearly shows three regions: a central object, the hot corona characterised by a temperature spike, and the 
accretion disk. The vertical lines indicate the limit of the convective mixing region (solid representing "deep" mixing and dotted for "shallow" mixing in the hot 
corona). See the online edition of the journal for a colour version of this figure. 
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Fig. 2. — Hydrogen deficient star surface abundance determinations taken fromlJefFery et alJ ROlT) compared to our nucleosynthesis yields computed in tlie 
merger of a 0.4 M© and a 0.8 M© wliite dwarfs. Tlie axis labels, [X] and [Fe], correspond to the logarithmic abundances relative to solar for individual elements 
and for iron, respectively. Solid lines are our calculations assuming "deep" convective mixing, while dashed lines are for "shallow" mixing (see text). The lines 
are added to guide the eye since our calculations include only two metallicities: Z = 0.015 and Z = 10" . The asterisk symbol corresponds to solar abundances, 
while the diagonal lines intersecting solar values are the abundances expected if the solar values were scaled with metalicity. See the online edition of the journal 
for a colour version of this figure. 



